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ABSTRACT 

Gravitationally unstable disks can fragment and form bound objects provided that their coohng time 
is short. In protoplanetary disks radiative coohng is Ukely to be too slow to permit formation of planets 
by fragmentation within several tens of AU from the star. Recently, convection has been suggested 
as a faster means of heat loss from the disk but here we demonstrate that it is only marginally more 
efficient than the radiative cooling. The crucial factor is the rate at which energy can be radiated from 
the disk photosphere, which is robustly limited from above in the convective case by the adiabatic 
temperature gradient (given a certain midplane temperature). Thus, although vigorous convection 
is definitely possible in disks, the inefficiency of radiative loss from the photosphere may create a 
bottleneck limiting the ability of the disk to form self-gravitating objects. Based on this argument 
we derive a set of analytical constraints which diagnose the susceptibility of an unstable disk to 
fragmentation and show that the formation of giant planets by fragmentation of protoplanetary disks 
is unlikely to occur at distances of tens of AU. At the same time these constraints do not preclude the 
possibility of fragmentation and star formation in accretion disks around supcrmassive black holes. 

Subject headings: planets and satellites: formation — solar system: formation — planetary systems: 
protoplanetary disks 



1. INTRODUCTION. 

It has been known for a long time (Safronov 1960; 
Toomre 1964) that massive differentially rotating disks 
are prone to gravitational instability (hereafter GI) which 
sets in when the destabilizing effect of the self-gravity in 
the disk exceeds the combined restoring action of the 
pressure and Coriolis forces. Instability condition is con- 
veniently expressed in terms of the so-called Toomre Q- 
parameter as 



Q 



< 



ttGS 



Qo, 



(1) 



where Vl, Cs, and S are the local angular frequency (disk 
is assumed to be in Keplerian rotation) , sound speed, and 
surface density of the disk, and Qo ~ 0.75 — 1.5 (Kim, 
Ostriker, & Stone 2002; Boss 2002) is a threshold value 
of Q at which GI sets in. 

Nonlinear outcome of the GI depends on the thermo- 
dynamics of the disk, namely on its ability to cool rapidly 
(Gammie 2001; Johnson & Gammie 2003). Rapid cooling 
keeps pressure forces subdominant compared to the disk 
self-gravity in the non-linear stage of GI thus promoting 
the collapse of unstable parts of the disk and leading to 
the formation of bound objects. When the cooling is slow 
fragmentation is avoided and the disk settles to a quasi- 
steady regime of gravitoturbulence (Gammie 2001) char- 
acterized by dramatic perturbations of the disk surface 
density and fluid velocity (Durisen et al. 2006). Gammie 
(2001), Rice et al. (2003), and Rice, Lodato, & Armitage 
(2005) have demonstrated using simulations with fixed 
cooling time tcooi that fragmentation requires 



(2) 



where ^ w 3 — 12 depending on the adiabatic index 7 of 
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the disk material, although in some situations C can be 
much higher (Johnson & Gammie 2003). 

Fragmentation resulting from GI is likely an impor- 
tant ingredient of the dynamics of accretion disks in the 
central regions of galaxies, around supcrmassive black 
holes (SMBHs; Paczynski 1978; Kolykhalov & Sunyaev 
1980; lUarionov & Romanova 1988; Levin 2003; Good- 
man 2003;Goodman & Tan 2004; Thompson, Quataert, 
& Murray 2005), including the past disk around the black 
hole in the center of our Galaxy (Levin & Beloborodov 
2003; Nayakshin 2005). On smaller scales, fragmenta- 
tion of gravitationally unstable protoplanetary disks re- 
sulting in direct formation of gas giant planets has been 
put forward (Cameron 1978; Boss 1998) as an alterna- 
tive to the conventional core instability model (Perri & 
Cameron 1974; Harris 1978; Mizuno 1980). Some simula- 
tions incorporating radiation transfer (Boss 2002; Mayer 
2006) seem to demonstrate that planet formation by GI 
is possible. However, these results are in drastic conflict 
with the long radiative cooling times of optically thick 
protoplanetary disks (Rafikov 2005; hereafter R05) and 
the recent numerical results by another group (Cai et al. 
2006; Boley et al. 2006). 

To reconcile inefficient radiative cooling with fragmen- 
tation exhibited in their simulations Boss (2004) and 
Mayer et al. (2006) have come forward with a sugges- 
tion that gravitationally unstable protoplanetary disks 
can be efficiently cooled by convection so that fragmen- 
tation becomes possible. The goal of this paper is to 
critically revise this possibility by setting in S^l a lower 
limit on the cooling timescalc of the disk in which en- 
ergy is transported away from the midplane by convec- 
tive motions. We then use this result in 21 to constrain 
the properties of the disks in which GI and convective 
cooling lead to fragmentation. 



2. CONDITIONS FOR CONVECTION. 
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The question of whether energy in the disk is trans- 
ported by convection or radiation depends on the behav- 
ior of disk opacity. Bell & Lin (1994) have demonstrated 
that in some temperature intervals opacity k{P, T) can 
be well represented as 



(3) 



{P and T are the gas pressure and temperature) which 
we assume to hold throughout the paper. 

In this section we focus on optically thick disks, for 
which r = J_^pKdz 3> 1. Following the approach of 

Lin & Papaloizou (1980) we demonstrate in Appendix 
El that some optically thick disks are guaranteed to be 
convective whenever 
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Condition directly follows from some assumptions 
about the spatial distribution of the energy dissipa- 
tion rate within the disk (see Appendix ^ and the 
Schwarzschild criterion for convection 



V > Vad, V = 



dliiT 
dlnP' 



(6) 



where V is a temperature gradient. 

Using condition Q one can infer that in the low- 
temperature regime disk is convective for T < 150 K, 
since then k is dominated by ice grains (a = 0, /3 « 2; 
Bell & Lin 1994) and V^d = 2/7 while Vq « 1/2 
(7 = 7/5) in this case. At higher temperatures, when ice 
grains sublimate (but below the sublimation point of the 
silicate and metal grains) a = 0, « 1/2 and Vq ~ Vad, 
so that the disk is marginally unstable for convection. At 
the same time, submillimeter observations of protoplan- 
etary disks demonstrate that dust opacity at low tem- 
peratures is best described by a = 0,/3 ~ 1 (Kitamura 
et al. 2002), most likely as a result of dust grain growth. 
In this case Vq = 1/3 and disk is convective. 

At higher temperatures, when all dust grains subli- 
mate, opacity is pressure-dependent (a ^ 0). In partic- 
ular, a w 2/3 and (3 « 7/3 in the regime of molecular 
opacity (T « 1.5 — 5 x lO'^ K), so that Vq ~ 1 and disk 
is again convective. When using criterion Q at even 
higher temperatures one has to bear in mind possible 
importance of the radiation pressure which we neglect in 
our analysis. 

3. UPPER LIMIT ON THE CONVECTIVE COOLING RATE. 

We first consider optically thick disks in which dense 
midplane regions cool in a two-stage fashion. Their heat 
first needs to be transported somehow to the disk pho- 
tosphere (located at r 1) and then lost to space by 
radiation from the photosphere, irrespective of whether 
transport of energy inside the disk is done by radiation 
or convection. 

In a convective disk the first leg of the cooling pro- 
cess can be rather fast but not instantaneous. Indeed, 
convection is driven by the buoyancy of the hot gas, but 
buoyancy cannot supply gas with the vertical acceler- 
ation higher than vertical gravitational acceleration gz- 



Because of that, when the parcel of hot gas reaches about 
a scale height above the midplane (which is roughly the 
location of the photosphere) its vertical speed would be 
at most ~ Cg and it would take time at least ^ 
to get there. Thus, the heat transport inside the disk 
cannot occur in less than the dynamical timescale of the 
disk. If this were the end of story, criterion ||2Jl would 
suggest that convective energy transport is marginally 
compatible with fragmentation. 

However, the need to eventually lose the transported 
energy by radiation from the photosphere changes this 
conclusion, since the cooling time is in fact the sum of 
the interior heat transport timescale (which can be as 
short as fi"^) and the timescale of radiative losses from 
the disk surface. For a given midplane temperature 
the fastest cooling occurs for the highest possible photo- 
spheric temperature Tph, since the energy loss from the 
photosphere per unit area is uT^f^, assuming blackbody 
emission. For a fixed higher Tph implies smaller tem- 
perature variation for a fixed change in pressure, or, in 
other words, lower value of V [see equation ©]. Thus, 
disk cools most efficiently when V assumes its lowest 
possible value. 

Whenever the radiative temperature gradient Vr de- 
fined as 

^^-16^ 

(here ct is a Stephan-Boltzmann constant and F is the 
upward energy flux) satisfles the condition Vr < Vad, 
disk is convectively stable and V = V,- < Vad- In this 
case energy is transported from the midplane of the disk 
to its photosphere radiatively, and equations ^ and l(2Jl 
directly yield constraints on the properties of disks ca- 
pable of forming bound objects by the GI (R05). As V^ 
becomes larger than Vad convection sets in and the tem- 
perature gradient obeys Vad < V < Vr. From this one 
can immediately see that the smallest possible value of 
V (leading to highest Tph) in a convective disk is Vad- 
Situation when V « Vad is possible only in the case of 
very efficient convection when the radiative flux is small 
compared to the convective flux, which typically requires 
high gas density. Whether this is the case depends on the 
situation at hand but using V = Vad allows us to obtain 
a robust upper limit on the cooling rate of the disk. 

Disk with V = Vad is vertically isentropic and we as- 
sume this to hold all the way to the photosphere (this is 
verified in Appendix . Then 



Tph 

Trn. 



Pph 



(8) 



where Pph and Pm are the gas pressure at the photo- 
sphere and at the midplane, respectively. For the pur- 
poses of this calculation it is enough to assume that ver- 
tical acceleration gz is roughly constant in which case 

Pph ^ ^i^PmiTm) 
Pm ^{Pph^Pph) 

Then one finds from @, lO, and © that 

Tph ^ „^/4 4Vad 

T 

m. 



(9) 



l + a + [iVad ^^^^ 
Thus, the fastest possible rate L of energy loss per unit 
surface area of an optically thick convective disk is i ~ 



3 



(jT^t^'^. Cassen (1993) used similar line of arguments to 
calculate Tph for the disk with constant opacity (which 
is, in fact, convectively stable). Equation (fTUIl reduces to 
his result i] = 4(7 — l)/7 if we set a = /3 = 0. 

In the optically thin case {t < 1) disk loses energy 
to space by volumetric radiative cooling irrespective of 
whether it is convective or radiative. Then, according to 
the Kirchhoff's law, L w <jT^t. In general, for arbitrary 
T, one can interpolate between the optically thick and 
thin regimes so that the cooling rate per unit surface 
area is 

fcir)^Xr^ + -, (11) 

where x, (/) ~ 1 are constants. In Appendix^we present 
detailed calculation of x and show that its value does not 
deviate strongly from unity. 

Cooling time is a ratio of the total thermal energy 
content per unit surface area of the disk Eth ~ Sc^ 
[cs = {kTm/ fJ^y^^ is the isothermal sound speed corre- 
sponding to the midplane temperature T^, the same as 
used in here fc is a Boltzmann constant, /i is the mean 
molecular weight] to the cooling rate L. Using 1)11(1 one 
finds the shortest possible cooling time in a convective 
disk: 



icool — 



■fc{r) 



(12) 



; 2 X 10* yrs 



T 



/c(r) 



103g cm-2 VlOO Ky 103 ' 

where the numerical estimate is for the gas of Solar com- 
position and one has keep in mind that /c(t) ~ 10^ 1 
for the adopted fiducial values of E and T. 

When the disk is convectively stable energy is trans- 
ported by radiation, however tcooi is still given by the ex- 
pression ((T^ but with fc replaced with /^(t) « t -|- 
(R05). In the optically thin case cooling is the same for 
both modes of heat transfer while in the optically thick 
regime one finds 



/c(r) 



(4-/3)(Vo-V,d) 



(13) 



Thus, when the condition Q is satisfied and the disk is 
convective, fastest possible convective cooling is more ef- 
fective than the radiative cooling would be for the same 
disk parameters (S, Tm,f2). However, this improvement 
is not very large for two reasons. First, in real disks only 
the condition Vad < V < weaker than V = Vod, 
must be fulfilled, so that the real cooling rate is some- 
where in between the radiative and the maximum possi- 
ble convective regimes. Thus, tcooi is likely to be longer 
than equation (|12|l predicts. Second, even the short- 
est cooling time H12f) is not so different from radiative 
tcooi- for example, ^ = 3/11 for a — 0, [3 — 2, while for 
a = 0, = 1 (as supported by observations of T Tauri 
disks) one finds ^ = 1/9. Thus, the difference between 
the radiative and convective tcooi is a very weak function 
of r: for r = 10^ the ratio of cooling times in two cases 
are fc/ fr ~ 7 and « 2 correspondingly. 

4. IMPLICATIONS FOR DISK FRAGMENTATION. 

Combination of gravitational instability condition ^ 
with the cooling requirement Q allows one to set a ro- 
bust lower limit on S and T of the disk region which can 



fragment as a result of GI, which has been previously 
done in R05 for the case of radiative heat transfer. Re- 
peating his arguments for the case of maximally efficient 
convective cooling [i.e. using eq. H12|l ] we find that self- 
gravitating objects can form only in disks satisfying the 
following properties: 



where 



I] > I]„ii 
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[fc{r)] 



1/5 



T»/[/c(T)]^/^ 
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(15) 



1 



-,1/5 



7/10 



6.6 X lO^g cm-a^S^/^" (Q[SA^C)"^' ) (1,6) 
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(17) 



Here oau is the semimajor axis a in units of AU, 
p, = fi/mn is the mean molecular weight relative to the 
atomic hydrogen mass mn, and Mi, is the mass of the 
central star {Mq is the Solar mass) . Needless to say, the 
condition H15|l has to be satisfied simultaneously with ^ 
which additionally limits T from above. 

One can clearly see that the only difference between the 
radiative case studied in R05 and convective case consid- 
ered here is in the explicit form of the function fcij), 
which changes to fr{T) in the radiative case. Bearing 
in mind that /c(t) is similar to fr{T) in that it reaches 
the minimum value /c ~ 1 at r '--^ 1 we conclude that the 
lowest possible {infimum) values of S and T given by H16|l 
and H17|l are the same as in the radiative case. Thus, all 
the constraints set in R05 on the basis of Si„ t and T, 



persist in the convective case as well. In particular, it 
follows from ^ and (0 that Tpu > T,„f[fc{T)f^^°, so 
that the photospheric temperature of fragmenting disk 
Tph is limited from below by at least Tm/. 

5. APPLICATION TO PROTOPLANETARY DISKS AND OTHER 
ENVIRONMENTS. 

For molecular gas of solar composition and C = 
10, Qo = 1 one finds Ei„/ « 2 x 10^ g cm^^ and 
Tinf ~ 1400 K at 1 AU. These very extreme conditions 
[E exceeds the corresponding value in the minimum mass 
Solar nebula (Hayashi 1981) by ^ 10^!] are clearly in- 
compatible with the observations of T Tauri disks which 
precludes the possibility of giant planet formation by GI 
within several AU from the parent star. At larger dis- 
tances, however, the constraints set by Ein/ and Tmf re- 
lax and become more acceptable from the observational 
point of view. At 10 AU, for example, Ej„/ « 1.6 x 10^ 
g cm^2 and T^f « 100 K which is at (or above) the 
uppermost end of the observed distributions of E and T. 

However, in such environment /c(t) « r ~ 10"^ ^ 1 
(for K « 1 cm^ g~^) so that equations (|14|I - H15|) set much 
more severe constraints on the disk properties. Careful 
inspection of 1(1411 - 1(15(1 demonstrates that planet forma- 
tion by GI at 10 AU is only possible if T > 1500 K so that 
dust has sublimated lowering the opacity — otherwise r 
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would be so high that according to H15() T cannot self- 
consistently remain below the sublimation point. Thus, 
planet formation by GI even at 10 AU from the cen- 
tral star requires disk properties which are clearly incon- 
sistent with the current observations of protoplanetary 
disks (Kitamura et al. 2002) . Only beyond 20 AU do we 
find that < 1500 K and Emin < 2.4 x 10^ g cm"^ 

so that silicates condense. At 45 AU T < 150 K and ice 
grains condense which increases opacity by a factor of 
several. In general, we find that the disk properties lim- 
ited by H14l) - (|15|l become more or less compatible with 
the observations of T Tauri disks only at a ~ 100 AU. 
Only this far from the star can planets possibly form by 
the GI, which contradicts the conclusion of Boss (2006). 

Previous discussion focused on the local disk properties 
necessary for planet formation by GI. It is however clear, 
by analogy with R05, that the global parameters of the 
fragmenting disk - its mass and luminosity - should also 
be rather extreme. We find that the disk mass would 
exceed several Mq if planets were to form by GI at 20 
AU, but it reduces to ~ 0.1 Mq at 100 AU. The masses 
of self-gravitating objects that would form as a result of 
fragmentation are also going to be large, likely in the 
brown dwarf regime (R05) , although this can be reliably 
clarified only by numerical simulations. 

Fragmentation of gravitationally unstable disks around 
SMBHs in galactic centers is also constrained by the con- 
ditions ((ni) - ((T3|l if these disks are convective. In this 
case, unlike the situation with the protoplanetary disks, 
no obvious conflict with the observations arises. Indeed, 
let's consider the possibility of star formation in a disk 
at a = 0.1 ps from the SMBH in the Galactic Center. 
For a black hole mass of 3 x 10^ Mq simultaneous frag- 
mentation and GI at 0.1 pc require Sm/ « 6 g cm~^ 
and Tinf « 4 K. The real temperature is in fact limited 
from below by the stellar irradiation at the level 40 K 
(Nayakshin 2005) and this automatically increases min- 
imum E to 25 g cm~^ to satisfy the Toomre crite- 
rion 1^. This translates into r « 8 and cooling time 
w 0.2n^^ which would lead to fragmentation. Charac- 
teristic disk mass tto^H, ~ 4 x 10'^ Mq is fully compatible 
with the current constraints on the total stellar mass in 
the two stellar disks around the central black hole in the 
Galactic Center (Nayakshin et al. 2006). 

6. DISCUSSION. 

There is good agreement between the constraints on E 
and T presented in this paper and R05 and the results 
of the three-dimensional grid-based simulations by Mejia 
(2004), Cai et al. (2006) and Boley et al. (2006). Using 
radiation transfer scheme based on the flux-limited diffu- 
sion these authors find that planet formation by GI does 
not occur at distances of ^ 10 AU from the star, in ac- 
cord with our findings. This outcome is clearly caused 
by the long radiative cooling timescales in dense and cold 
disk models that are run by these authors. In particular, 
Cai et al. (2006) find that lowering disk metallicity has 
the effect of decreasing the cooling timescale, which is 
obvious from equation (|12|l if one also notices that the 
low-temperature opacity n is proportional to metallicity 
and /c(t) oc in the optically thick regime (typical at 
distances of several tens of AU in massive protoplanetary 
disks). 

At the same time recent radiative hydrodynamical sim- 



ulations by Boss (2004) and Mayer et al. (2006) demon- 
strate efficient fragmentation of massive and cold proto- 
planetary disks and thus disagree with both the analyt- 
ical constraints presented in this paper and the simula- 
tions of Mejia (2004), Cai et al. (2006) and Boley et al. 
(2006). The arguments presented here do not support the 
idea of Boss (2004) and Mayer et al. (2006) that energy 
transport by convection can significantly speed up the 
disk cooling and lead to its fragmentation. While there 
are many factors that can cause this discrepancy (opac- 
ities, treatment of shock dissipation, numerical scheme 
employed, etc.) we believe that the major source of dis- 
agreement is in the treatment of the external radiative 
boundary condition by different groups. 

To illustrate this point we note that Boss (2004) , Mayer 
et al. (2006), and Boley et al. (2006) do observe^ efficient 
convection in their simulations. However, rapid vertical 
transport of energy from the midplane to the disk sur- 
face is a necessary but not sufficient condition for rapid 
cooling, since this energy must then be radiated away 
from the photosphere. Thus, the photosphere presents a 
bottleneck for the disk cooling such that the disk cannot 
cool faster than its photospheric temperature permits. 
For this reason it is very important that the external ra- 
diative boundary condition is treated with extreme care^ 
in simulations. 

Some of the vertical motions seen in simulations and 
attributed to convection may actually be due to the up- 
ward propagating waves launched by the overdensities 
caused by the GI (Boley et al. 2006). The dissipation 
of such waves may heat up the upper layers of the disk 
which has an effect of increasing the photospheric tem- 
perature. This, however, has nothing to do with the 
cooling of transiently collapsing overdensities in the disk 
as these waves do nothing to transport the thermal en- 
ergy away from the dense gas near the midplane. Thus, 
such wave motions cannot facilitate fragmentation of an 
unstable disk. 

Similarly, the external irradiation of the disk (by the 
central star or the surrounding envelope) which acts to 
increase the photospheric temperature does not acceler- 
ate disk cooling. In this case, although the disk surface 
loses more energy because of the higher Tph it also gains 
energy from the absorbed external radiation which keeps 
the net cooling rate unchanged. Cooling of overdensities 
forming near the midplane as a result of GI depends only 
on the temperature gradients that develop in the interior 
and thus is not affected by external heating. Besides, 
external illumination tends to stabilize temperature gra- 
dient and to suppress convection. 

Our final comment refers to the use of constant fac- 
tor C ~ 1 — 10 in the fragmentation condition Using 
simulations with nonlinear cooling rates devised to mimic 
the effect of realistic opacities on disk cooling Johnson & 
Gammie (2003) have demonstrated that C can be signif- 
icantly larger than ~ 10 just below the so-called opacity 
gaps — regions in the p—T space where opacity suddenly 
changes. In particular, they find that C, ^ 10** at T ~ 10^ 

^ Boley et al. (2006) find convection to be disrupted during the 
active phase of the GI. 

■* In particular, the proper treatment of the photospheric tran- 
sition between the optically thick and thin regions should be espe- 
cially challenging for the SPH simulations such as used by Mayer 
et al. (2006). 
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K, near the point of silicate dust sublimation where opac- 
ity drops by a factor of ^ 10'^. Such an increase of ( sig- 
nificantly relaxes our constraints l|16|) - (|17|l which seems 
to allow fragmentation in disks of lower E and T. How- 
ever, the less conservative (but still robust) constraint 
I|14(I - (I15|I eliminates this concern because (similar to C) 
K and T are also very large on the verge of the opacity gap 
(compared to their values at the bottom of the gap) so 
that /c ^ 1 in (|14|I - H15|I and this offsets the decrease of 
"Sinf and Tinf driven by large C. Besides, this issue only 
arises if T ~ 10'^ K which is atypical for a protoplanetary 
disk in T Tauri phase anyway. 

The arguments presented in this paper and R05 are 
analytic in nature, however, they are based solely on the 
numerically verified criteria and Thus, keeping 
in mind previous discussion, we conclude that to within 
factors of order unity our constraints H14|l - (|15|l on the 
properties of unstable disks capable to support the for- 
mation of self-gravitating objects should be valid in ap- 
plications to realistic disks in which energy is transported 
by convection. 

7. CONCLUSIONS. 

We have investigated the possibility of rapid disk cool- 
ing by convection in the context of fragmentation of grav- 
itationally unstable disks. We have shown that even 
the most extreme form of convective cooling (realized 



when V — Vad) produces cooling rates which are only 
marginally higher than in the case of radiative energy 
transport. The major reason for the inefficient cool- 
ing of dense and cold convective disks is in the bot- 
tleneck caused by the need to eventually radiate from 
the photosphere the energy that has been transported 
there from the midplane, irrespective of how the latter 
has been done. Armed with this knowledge we demon- 
strate that inefficient cooling precludes direct formation 
of giant planets by GI in protoplanetary disks anywhere 
within several tens of AU from the parent star (analo- 
gous conclusion has been reached in R05). At the same 
time, disk cooling (convective or radiative) is fast enough 
to allow fragmentation of gravitationally unstable disks 
around black holes in the galactic centers, supporting the 
idea that the efficient star formation is possible in such 
disks. 
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APPENDIX 

DERIVATION OF CONDITION (gj . 

To derive the condition ^ in the case of opacity given by equation Q we will follow the approach adopted by the 



Lin & Papaloizou (1980). Equation of hydrostatic equilibrium states that 

dP 

dz 



-P9z, 



(Al) 



where vertical acceleration = ^^z + g^^ is contributed both by the central object {VPz term) and the disk self-gravity 
[gl^ term). The latter contribution satisfies Poisson equation 



Radiative transport is described by 



dz 



IGcrT^ dT 



dz 



(A2) 



(A3) 



3kp dz 

where F is the vertical radiative flux and e is the volumetric energy dissipation rate. Equations IjAip and ljA3|l result 
in Vr given by equation (jT)), but they can also be combined in the following way: 

16 a d{T*-^^) _ F 

gz' 



Vo- 



(A4) 



3 Kod(pi+") 

We now note that if F/gz monotonically decreases as z increases (and, accordingly, as P decreases), then, integrating 
equation (jA4|) down from the photosphere we arrive at the following inequality: 

p 



^ 16 cr / 
Vo-- ' 

3 Kn V 



3 kq 



T 



-4-/3 



T. 



4-/3 



ph 



J gz 



< - P 



- P 



1 + Q 

ph 



(A5) 



since F/ g^ is assumed to be largest at P. This can be finally rewritten as 



Vc 



1 



Tph 



4-/3' 



< V, 



1 



Pph 
p 



1+a 



(A6) 



which reduces to Vo < at high depth in the disk, where T ^ Tp^ and P ^ Pp^. Thus, if (1) the condition is 
fulfilled and (2) F/g^ is a decreasing function of z, then Vr > Vad and convection has to operate in the bulk of the 
disk since, according to the Schwarzschild criterion ©, disk cannot be radiative. 
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We now find out under which circumstances F / monotonically decreases as z grows. One has 

z z 

§Tz{^)^J ^'^'^ '^''^^ + ^ / ^<''^P^'''^ <^')P(^)] dz', (A7) 



where equation (|A2p has been used. The first term in the right-hand side is clearly negative when a reasonable and 
weak assumption of e decreasing with increasing z is made. Sign of the second term, which quantifies the effect 
of the disk self-gravity on stability, depends on how rapidly e decreases with z. If, like in a conventional a-disk 
(Shakura & Sunyaev 1973) with constant a and sound speed, e{z) oc p{z)^ then the second term is identically zero 
and d{F/gz)/dz < 0. When e(z) decays with z faster than p{z) the second term in (jA7p is negative an d d{F /gz)/dz is 
again negative. However, if e{z)/ p{z) is an increasing function of z then the sign of the second term in l|A7p is positive 
and one cannot tell for sure how F/gz behaves as a function of z. In the latter case condition Q may or may not be 
sufficient to determine whether disk is convective. 

It may be possible that gravitationally unstable disks have e{z)/p{z) increasing with z. This can be caused, for 
example, by dissipation of the upward-propagating waves driven by the instability, which would tend to deposit their 
energy high up in the disk stabilizing its structure. We can then conclude that in many situations (certainly in non- 
self-gravitating disks) condition is sufficient for convection, while in some self-gravitating disks a more conservative 
criterion than Q may be required. 



COOLING OF AN ISENTROPIC DISK. 

Here we calculate the structure and the cooling time of a disk that has an isentropic interior (for which the condition 
V = Vad is satisfied) with equation of state P = Kp^ smoothly matching the outer radiative layer near the photosphere. 
We will separately consider structure of non-self-gravitating disks with = fi^z and self-gravitating disks with gz = gl^ 
satisfying equation (|A2p . 

Non-self-gravitating disk. 

In the non-self-gravitating isentropic disk equation (|A1|I can be solved using g^ = fi^z, resulting in the following 
disk structure: 



T{z) 
T 



P(z) 



7-1 



^2' 



2 ksTrr 



(Bl) 



where Tm, Pm, Pm are the midplane temperature, pressure, and density, while H is recognized as the height of the disk 
surface. Using I|B1|I and assuming disk to be optically thick (so that the isentropic part contains most of the mass) 
one can calculate the total surface density 



and the full optical depth 



oo i 

Y, = 2jpdz^ hpmH, /i = 2 y (1 - x^/^^-^Ux, 



oo 1 

r = 2j pKdz = hPrn^i^H, h = 2 a;')'3+("^+i)/(^-i)dx, 



(B2) 



(B3) 



where Hm = K{Pm, Tm) — KoP^T^ is the opacity at the midplane. The total thermal energy content per unit surface 
area is 

oc 1 



Eth — 



7-1 



Pdz = hPrnH, h = 



7-1 



{l-x'^Y'^-^dx, 



(B4) 



while the midplane pressure P„i — (^ad/'^)^'^H'^Pm [following from eq. ljBl|) and the ideal gas law] can be expressed 
through the value of g^ at the disk surface gz{H) = (l^H in the following way: 



Pm=hPrnHgziH), h^^- 



(B5) 



Self- gravitating disk. 

In a self-gravitating case one can solve the system (|Aip - (|A2(l with the equation of state P — Kp'' to find the following 
implicit relation between p and z: 

1 

z^HIq^ J x^-^{l~x^y^/'^dx, (B6) 

P/Pm. 



where 



1/2 



Jo, h= / x'<-^l-x'<)-^'^dx, 



(B7) 



K 

is the height of the disk surface. Using (|B6p we can derive expressions for S, t, Eth, Pm identical to equations l|B2|) - (|B5l 
but with constant coefficients /fe, A; = f , .., 4 which are now given by 



a7+(/5+l)(7-l)(l_^7)-l/2j^^ 



2Jn 



7-1 



2.2(7-l)(l_^7)-l/2^2.^ 



(B8) 

(B9) 

(BIO) 
(Bll) 



To derive (|Bf fp one should notice that gz{ii) = while P™ = A'p^Jj = 8ttGpI^H^{jIo)^^ , see equation (|B7(I . 

Outer radiative layer. 

Somewhere near the disk surface, at z « iJ, a transition from the convective interior to the radiative atmosphere 
occurs. We assume that this happens at the temperature Ttr and pressure Ptr- Inside the radiative layer 

(B12) 



Differentiating this equation with respect to z and using equation IjAf p with constant gz{H) (which is justified by the 
small thickness of the radiative layer compared to H, following from the assumption Tph <^ Tm) one finds that 



2(/3-4)/4 ^ X 



1/(4-/3) 



X 



T — Tph 

in the outer radiative zone, so that 

oo 

T[z)=jKpdz=^ (2(^-4)/4 + x) 

Z 

[here we have used equation (|Al|l ] and 



16Vog.(H) 



-+a 



4/(4-/3) 1 

^ 2 



V. = ^^^^ = Vo- 



X 



d\nP ^'2('3-4)/4 + x' 
Transition from radiative to convective energy transport occurs when = Vad, i-e. at 

_ 2(/3-4)/4V,, 
(Vo-V,,)' 

while the photosp here ( r = 2/ 3) lies at X ph — 1 — 2^^~^^^^, see equation (jB14|l . 
From equations (|B1|) . I|B13|I . and (|B16() one finds 



(B13) 



(B14) 



(B15) 



(B16) 



Ttr = 



Ptr — 



Tph ( Vo 



21/4 VVo- Vad 



1/(4-/3) 



16x2(/^-4)/4v„v^,g^(g) 



1/(1+q) 



Using these relations, equations HB3|) and HB5() . and the fact that Ttr/T„i — {Ptr / Pm)'^ one obtains 



T 

J. 



/ l6Varf/2 

V 3/4 



Vad /o(/3-4)/4 



Vo 



Vo - Vad 



V„d-Vo' 



l/(l+Q+/3V„d) 



(B17) 
(B18) 

(B19) 
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where was introduced in equation H10|l . Using equations ljB3|) . (|B4p and (|B19|I we finally derive cooling time in the 
optically thick case as 

Uooi = ^ = X^r\ X = (B20) 

In p a rticu lar, in the case of a non-self-gravitating disk one finds for a = 0,/3 = 2 (77 = 8/11) using equations 
(|B2l-(lB5l) that A = 1.37, x ^ 0.31, while for a = 0, /3 = 1 ij ] = 8/9) one gets A = 1.55, x = 0.19. When the 
disk is self-gravitating one obtains using equations (|B8p - (|Bll|l that A — 1.24, x = 0.46 for a = 0,/3 = 2, while 
A — 1.37, X = 0.3 for a = 0,/3 = 1. In realistic unstable Keplerian disk the self-gravity of the disk and the central 
star contribute roughly equally to so that one should expect A = 1.24 — 1.37, x = 0.31 — 0.46 for a = 0, /9 = 2 and 
A = 1.37 — 1.55, X = 0.19 — 0.3 for a = 0,/3 = 1. Note that the small value of x in some cases should not strongly 
affect our estimates of E^in and Tmin as these quantities depend on x rather weakly, see equations (|14|1 - H15() . 

REFERENCES 



Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 

1990, AJ, 99, 924 
Bell, K. R. & Lin, D. N. C. 1994, ApJ, 427, 987 
Boley, A. C, Mejia, A. C, Durisen, R. H., Cai, K., Pickett, M. K., 

& D'Alessio, P., 2006, astro- ph/0607112 
Boss, A. P. 1998, ApJ, 503, 923 ' 
Boss, A. P. 2002, ApJ, 576, 462 
Boss, A. P. 2004, ApJ, 610, 456 
Boss, A. P. 2006, ApJL, 637, L137 

Cai, K., Durisen, R. H., Michael, S., Boley, A. C, Meja, A. C, 

Pickett, M. K., & D'Alessio, P. 2006, ApJL, 636, L149 
Cameron, A. G. W. 1978, Moon and the Planets, 18, 5 
Cassen, P. 1993, in "Lunar and Planetary Inst., Twenty-fourth 

Lunar and Planetary Science Conference", Part 1: A-F, 261 
Durisen, R. H., Boss, A. P., Mayer, L., Nelson, A. F., Quinn, T., 

& Rice, W. K. M. 2006, astro-ph/060317S3^ 
Gammie, C. F. 2001, ApJ, 553, 174 
Goodman, J. 2003, MNRAS, 339, 937 
Goodman, J. & Tan, J. C. 2004, ApJ, 608, 108 
Harris, A. W. 1978, Lunar Planet. Sci., 9, 459 
Hayashi, C. 1981, Progr. Theor. Phys. Suppl., 70, 35 
Illarionov, A. F. & Romanova, M. M. 1988, Sov. Astr., 32, 148 
Johnson, B. M. & Gammie, C. F. 2003, ApJ, 597, 131 
Kim, W.-T., Ostriker, E. C, & Stone, J. M. 2003, ApJ, 595, 574 
Kitamura, Y., Momose, M., Yokogawa, S., Kawabe, R., Tamura, 

M., & Ida, S. 2002, ApJ, 581, 357 
Kolykhalov, P. I. & Sunyaev, R. A. 1980, SovAL, 6, 357 



Lin, D. N. C. & Papaloizou, J. 1980, MNRAS, 191, 37 
Levin, Yu. 2003, astro-ph/0307084 

Levin, Yu. & Beloborodov, A. M. 2003, ApJL, 590, L33 

Mayer, L., Lufkin, G., Quinn, T., & Wadsley, J. 2006, 

astro-ph/0606361 
Mayer, L., Qumn, T., Wadsley, J., & Stadel, J., 2002, Science, 298, 

1756 

Mizuno, H. 1980, Progr. Theor. Phys., 64, 544 
Nayakshin, S. 2005, astro-ph/0512255 

Nayakshin, S., Dehnen, W., Cuadra, J., & Genzel, R. 2006, 

MNRAS, 366, 1410 
Paczynski, B. 1978, AcA, 28, 91 

Perri, F. & Cameron, A. G. W. 1974, Icarus, 22, 416 
Pickett, B., Cassen, P., Durisen, R. H., & Link, R. 1998, ApJ, 504, 
468 

Rafikov, R. R. 2005, ApJL, 621, L69 (R05) 

Rice, W. K. M., Armitage, P. J., Bate, M. R., & Bonnell, I. A. 

2003, MNRAS, 339, 1025 
Rice, W. K. M., Lodato, G., & Armitage, P. J. 2005, MNRAS, 364, 

L56 

Safronov, V. S. 1960, Ann. d'Astrophysique, 23, 979 
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
Thompson, T. A., Quataert, E., & Murray, N. 2005, ApJ, 630, 167 
Toomre, A. 1964, ApJ, 139, 1217 



